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We study light propagation in nanoscale periodic structures composed of dielectric and metal in the visible 
range. We demonstrate that diffraction can be tailored both in magnitude and in sign by varying the geometric 
features of the waveguides. Diffraction management on a subwavelength scale is demonstrated by numerical 
solution of Maxwell equations in frequency domain. © 2008 Optical Society of America 
OCIS codes: 050.6624, 250.5403, 130.2790, 160.3900 



The miniaturization of photonic devices for confining 
and guiding electromagnetic energy down to nanometer 
scale is one of the biggest challenges for the informa- 
tion technology industries. 1 In the last years, photonic 
crystals technology allowed to gain one order of magni- 
tude factor in the miniaturization of components such 
as waveguides and couplers with respect to conventional 
(i.e. based on total internal reflection) optics. However 
when the size of a conventional optical circuit is reduced 
to the nanoscale, the propagation of light is limited by 
diffraction. One way to overcome this limit is through 
surface plasmon polaritons, 2 which are evanescent waves 
trapped at the interface between a medium with positive 
dielectric constant and one with negative dielectric con- 
stant, such as metals in the visible range. Even though 
this phenomenon has been known for a long time, in the 
last years there is a renewed interest in this field, mainly 
motivated by the will to merge integrated electronic cir- 
cuits to photonic devices. 3 

In this Letter we study the propagation of light in 
nanostructured metal-dielectric waveguide arrays (plas- 
monic arrays). As well known, an array of evanes- 
cently coupled single-mode waveguides exhibits power 
exchanges among the waveguides leading to discrete 
diffraction. 4 ' 5 In plasmonic arrays we find peculiar 
diffractive phenomena and we show the possibility of 
diffraction management on a subwavelength scale. 

The basic building block of a uniform waveguide ar- 
ray is the directional coupler; in the following, for the 
sake of clarity, we will consider only two dimensional 
cases: we have translational invar iance along the z axis, 
with y being the propagation direction and x the other 
transverse coordinate. In its simplest form, a directional 
coupler consists of two identical parallel waveguides in 
close proximity; as well known, the power exchange be- 
tween the two waveguides can be described by ordinary 
differential equations coupling the modal field amplitude 
^1,2(2/) of waveguides 1 and 2: 
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Fig. 1. (Color online) a) Dispersion relation of the funda- 
mental (odd) and second order (even) mode of the cou- 
pler DC1. b) Fundamental (continuous red) and second 
order (dash-dot blue) mode @ 600nm. 



where (3 is the propagation constant of the waveguide 
mode and C is the coupling coefficient, whose expression 
can be obtained in the framework of Coupled Mode The- 
ory. Limiting our attention to the TM case (i.e the non 
vanishing field components are H Zl E Xl E y ) we have an 
expression valid both for y- homogeneous and ^/-periodic 
(photonic crystal) waveguides: 6 ' 7 
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where e(x,y) is the permittivity of the directional cou- 
pler, £\, 2 {x,y) is the permittivity of the isolated single 
waveguide 1 or 2, e xli2 (x,y), e ylj2 (x,y) and h zli2 (x,y) 
are the electric and magnetic field of the mode (or the 
Bloch mode in case of periodic waveguide) of waveguide 
1 or 2, L y is the unit cell length for a photonic crystal 
structure and an arbitrary length in the case of transla- 
tional invar iance along y. 

In contrast with conventional optical devices based on 
total internal reflection where C must be always positive, 
in the following we introduce two different plasmonic di- 
rectional couplers having opposite sign of the coupling 
coefficient. 

Let us first consider a system composed of alternate 
layers of metal (for example Silver) and dielectric (air). 
We describe the optical properties of the metal using a 
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Drude free-electron model e(oj) = 1 — ujp/[uj(u — iu; r )] : 
where uj p is the plasma frequency and uo T is the colli- 
sion frequency. As an example we consider a structure 
composed of 30nm of Silver and 120nm of air. We ne- 
glect absorption in metal (lj t = 0) since it does not af- 
fect strongly the dispersive properties of the propagating 
modes. In Fig. [I] we show the dispersion relation and the 
profiles of the two modes supported by this directional 
coupler (to be called DC1 in the rest of the paper). In 
stark contrast with conventional waveguides, the funda- 
mental mode of this structure is odd. 8 Moreover we can 
see that the fundamental mode has one node whereas the 
second one has no. This feature seems to struggle against 
the well known Sturm oscillation theorem, 9 which states 
that the Nth order mode has exactly N — 1 nodes (N = 
1,2,...); however it has to be remembered that the the- 
orem holds only if the dielectric constant is everywhere 
positive. Note also that the reversal of modes parity 
implies a negative value of the coupling coefficient 10, 11 
C — (/3 even — /3 dd)/2- For example at a wavelength of 
600 nm Cab = (Peven - Podd)/2 = -6.59 • lO 5 ™" 1 and 
from Eq. (|f C C mt = -5.75- lO 5 ™" 1 . The difference be- 
tween the two values (around 13%) is due to the strong 
coupling between the waveguides. 

Another well known plasmonic guiding structure is 
the metal nanoparticle array, 12 where the energy trans- 
port is supplied by electromagnetic resonant coupling be- 
tween metal particles arranged in a linear chain. Double 
nanoparticle chains, where the electromagnetic energy 
is confined between two linear chains, offer a more flex- 
ible structure since the propagation is less determined 
by resonances, 13 allowing for a larger bandwidth for the 
guided modes. As an example we considered double chain 
waveguides composed of Silver nanoparticles with a ra- 
dius of 50 nm immersed in air with a longitudinal sep- 
aration of 110 nm and a distance between the chains 
of 150 nm. We neglected losses in the metal since it was 
found from previous band diagram calculations of metal- 
lic photonic crystals that this is reasonable for realistic 
absorption coefficients. 14 In Fig. ||]we show the disper- 
sion relation and the mode profiles of two coupled waveg- 
uides composed of three nanoparticle chains (to be called 
directional coupler DC2 in the rest of the paper). Since 
we are interested in guided modes we consider only the 
region in the (k,u) space that lies below the light line 
co = ck. This coupler supports several modes, however in 
the working range [450nm-750nm] (light shaded region 
in Fig. ) only two modes are reasonably below the 
light line. As opposed to the previous directional cou- 
pler, the fundamental mode here is even and the second 
order mode is odd, implying positive coupling constant. 
For example at 600nm we obtain Ca/3 = 1.29 • 10 6 m _1 
and Ccmt = 1-15- lO 6 ™" 1 . 

Let us now consider the uniform arrays Al (A2) ob- 
tained using directional couplers DC1 (DC2) as basic 
building blocks. The set of equations modelling the evo- 
lution of the modal field amplitude A n in each waveguide 




k-a/(2 7t) 



Fig. 2. (Color online) a) Dispersion relation of the cou- 
pler DC2. Thick curves, fundamental (even) and second 
order (odd) mode. Dashed curves, high order modes. 
Thin horizontal line denotes a wavelength of 600 nm. 
Dark shaded region, light cone; light shaded region, op- 
erating bandwidth 450nm-750nm. b) Fundamental and 
second order mode @ 600nm. 

of an uniform array reads as: 4 ' 5 
dA 

i-T- + P A n + C(A n+1 + A n _x) = 0. (3) 
dy 

The diffraction coefficient 15 ' 16 of the waveguide array 
can be derived from Eqs. 

D = -2Cd 2 cos(k x d), (4) 

where d is the spacing between the centers of adjacent 
waveguides and then k x d is the imposed input field phase 
shift between them. 

From the above reported analysis we expect, at nor- 
mal incidence (k x = 0), the waveguide array Al to be an 
anomalous-diffraction array (D > 0), whereas waveguide 
array A2 is expected to behave as a normal-diffraction 
array (D < 0). The field evolution along the waveguide 
arrays has been simulated without any approximation 
by solving Maxwell's equations through a frequency- 
domain finite-element method, using arrays composed 
by 17 waveguides. The central waveguide of the arrays 
is excited with a Gaussian field, which spreads during 
propagation and generates the typical diffraction pattern 
observed also in conventional waveguide arrays: two out- 
ermost wings and a few less intense peaks in the central 
waveguides. The same qualitative behavior is observed 
for both arrays since the intensity evolution is not in- 
fluenced by the diffraction sign. On the other hand the 
phase front curvature of the propagating field depends 
on the diffraction sign. Therefore, if we alternate arrays 
characterized by normal and anomalous diffraction, the 
input field shape can be periodically recovered. 

In Fig. [3] we report the finite element simulation of 
the diffraction-managed device. The first array section 
(Al) is Li=3900nm long, the second (A2) L 2 =1600nm. 
We excited only the central waveguide with a normal- 
incident 120nm FHWM intensity Gaussian field at 
600nm. The input excitation spreads during propaga- 
tion in the first section, whereas it exhibits an opposite 
behavior in the second one, and the initial field distribu- 
tion is recovered at the output end of the device. The 
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Fig. 3. (Color online) Time average power flow in y di- 
rection (normalized to the maximum) in the diffraction 
managed device @ 600nm. 



reported phenomena are a clear signature of the inver- 
sion of the sign of the diffraction coefficient, in perfect 
agreement with our theoretical analysis. It is worth not- 
ing that all the dynamics takes place in a 2/xra x bfim 
device. 

As far as losses are concerned, we verified that the 
propagation in the diffraction managed devices is not in- 
fluenced at all by including a lossy model for the metal 
(Drude model gives s(600nm) « — 15 — 0.37i): the evolu- 
tion showed in Fig. [3] is indistinguishable from the ana- 
logue calculated with a lossy metal. Moreover the decay 
length of the fundamental mode of the waveguides com- 
posing Al is, for example, L£>(600nm) = [2Im(/?)] _1 = 
17/im, much longer than the device length, indicating 
that all the relevant dynamics can take place without 
being suppressed by absorption. 

We now turn our attention to the bandwidth of the 
reported phenomena. By observing the band structure 
of the metal nanoparticle coupler (Fig. [2] a)), we ex- 
pect the bandwidth of this device to be the interval 
[450nm,750nm], i.e all the visible range from blue to red. 
We studied the dynamics of diffraction management in 
this range and we found that, even if the two guiding 
structures have very different dispersive properties, the 
diffraction compensation is reasonably good for all the 
wavelengths range. Figures j4^)-b) show the time aver- 
age power flow evolution at the edges of the operating 
frequencies interval. At short wavelengths (Fig. [4^)) the 
magnitude of the dielectric constant of metal is relatively 
low (« —8), the fundamental mode of the waveguides is 
poorly confined and the coupling is strong, leading to 
large diffraction in both arrays. Whereas at long wave- 
lengths (Fig. (4|d)) the large magnitude of the dielectric 
constant of metal (« —24) leads to strong confinement 
and small diffraction. The average value of diffr active 
spreading however remains low, considering the huge 
band we are looking at. Figure [4] c) displays the dis- 
persion of the diffraction parameter (coupling coefficient 
times propagation length C • L) for the arrays Al and 
A2 and for the entire device. 

In conclusion, diffraction properties of uniform arrays 
of plasmonic waveguides have been studied. Starting 
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Fig. 4. (Color online) Time average power flow in y di- 
rection (normalized to the maximum) in the diffraction 
managed device at different wavelengths : a) 450nm, 
b) 750nm. c) Diffraction coefficient times propagation 
length (C • L) for array Al (red), array A2 (blue) and 
for the entire device (C\ • L\ + C2 • £2) (black). 

from the analysis of coupling between adjacent waveg- 
uides we have demonstrated that diffraction can be con- 
trolled both in amplitude and sign with normal incidence 
input excitations. Diffraction management in an ultra- 
compact device composed of alternated sections with op- 
posite diffraction sign has been demonstrated on a wave- 
length interval covering the visible range from blue to 
red. 
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